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lO , The spin- 1/2 XXZ Heisenberg magnet is considered for the case of the anisotropy 

parameter tending to infinity (so-called, Ising limit). A thermal correlation function 
. of the ferromagnetic string is calculated over the ground state. The approach to 

[ the calculation of the correlation functions is based on the observation that the 

' wave function in the limit considered is expressed in terms of the symmetric Schur 

functions. It is demonstrated that at low temperatures the amplitude of the asymp- 
totical expression of the present correlation function is proportional to the squared 
numbers of strict boxed plane partitions. 
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1 Introduction 



The Ising model was one of the favorite models of A. N. Vassiliev. His studies were 
concentrated on the scaling properties of this model [HE]- It is remarkable that in a 
special limit of the anisotropy parameter tending to infinity, the Hamiltonian of the spin- 
1/2 Heisenberg XXZ chain is transformed into the Hamiltonian of one-dimensional Ising 
model [3J. This limit is usually called as the Ising limit of the XX Z chain [4J. In our paper 
we shall investigate the Ising limit of the XXZ magnet using an effective Hamiltonian 
which is commutative with the Ising Hamiltonian and has a complete system of the eigen- 
functions which is common with that of the Ising Hamiltonian. This effective Hamiltonian 
is an appropriate object for studying of the Ising limit since it is related to the transfer- 
matrix of the four-vertex model, which arises in the limit of infinite anisotropy from the 
transfer-matrix of the six-vertex model. In their turn, the trace identities relate the XXZ 
Hamiltonian and the transfer-matrix of the six-vertex model. 

The problem of calculation of the correlation functions of XXZ Heisenberg model 
in the framework of the Quantum Inverse Scattering Method |5J has required serious 
efforts [SHU]. We shall show that calculation of the correlation functions in the Ising 
limit is closely related to the theory of the symmetric functions [12j, and thus is closely 
connected with the problems of enumerative combinatorics, such as the random walks of 
vicious walkers |13ffI7] and enumeration of boxed plane partitions [T8|IT9]. In this paper 
we shall calculate the thermal correlation functions of the ferromagnetic string over the 
ground state using the approach developed in |20H24| . 

The XXZ spin- 1/2 Heisenberg magnet is defined on the one-dimensional lattice con- 
sisting of M + 1 sites labeled by elements of the set Ai = {0 < k < M, k G Z}, 
M + 1 = (mod 2). The corresponding spin Hamiltonian is defined as follows: 

1 M A 
#xxz = - g $^K+i4 + a t+i a k + ~2 (°Ui°l -l) + h°t)> (!) 

fc=0 

where the parameter A £ 1R is the anisotropy of the model, and h is the external magnetic 
field. The local spin operators = |(crf ±zcr^) and cr|, dependent on the lattice argument 
k G Ai, are defined standardly (see [6j[7]), and their commutation rules are given by the 
relations: [o'^,cr/] = S^af, t !; ^] = ±2 5^ cr^. The state-space is spanned over of 
the state- vectors (§)flo > where s implies either j" or I. The spin "up" and "down" 
states (|t) and \\), respectively) provide a natural basis in the linear space C 2 so that 

(J), m?). 

The periodic boundary conditions &t+(h[+i) = a t are imposed. The Hamiltonian (JTJ 
commutes with the operator of third component of the total spin S z : 

~ 1 M 

[H XXZ ,S Z ] = 0, S z =-J2*k- (2) 

fc=0 

To represent iV-particle state- vectors, | ^^(ui, . . . , un)), let the sites with spin "down" 
states be labeled by the coordinates fa forming a strict partition /x = /i 2 , • • • , A*jv), 
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where M > > pi^ > ■ ■ ■ > A*jv > 0. There is a correspondence between each partition 
and an appropriate sequence of zeros and unities of the form: je^ = efc(/x)} fcgyV) , where 
= 5fc i(Un , 1 < n < N . The Hamiltonian ([T]) is diagonalised via the following ansatz: 

M 

| X) xf Z (u) IlK) efc 1^) , (3) 

where the sum is taken over strict partitions /z, and = (^) n=0 |t)n is the fully 
polarized state with all spins "up". It is proposed to use bold-faced letters as short-hand 
notations for appropriate iV-tuples of numbers: for instance, u instead of (u\, . . . ,Un), 
etc. Therefore, the wave function X^ XZ ( U ) i n © is of the form: 

x? z (u) = E ^( u ) OlT ■ ■ ■ < N > ( 4 ) 

S P1 , P2 ,...,p N 

„ l-2Au=+ u 2 u 2 
As u = 1 2 P > 2 Pl Pl , (5) 

X<j<i<N Pj 



where summation goes over all elements of the symmetric group 

1, 2, N~ 
\pi, P2, PN; 



£>pi,P2,—,pn — 



The state- vectors ([3]) are the eigen-states of the Hamiltonian ([Q), if and only if the variables 
Ui satisfy the Bethe equations 

, 2(M+i) , 1 s N -i t~t 1 - 2Am 2 + m 2 m| 
' =( } M 1 - 2A«2 + « ' 1 - l ~ N - (6) 

The corresponding eigen-energies are given by 

1 - 

E nH = - J E( M ' + ^ 2 - 2A - 2/i) . (7) 

8=1 

The limit of zero anisotropy, A — > 0, transforms ([1]) to the Hamiltonian of the XX 
model on the periodic chain 

1 M ~ ~ 
H X x = -g E + ^ ' = ^fe+i^fe + ^fc+i * • ( 8 ) 

fc=0 

This limit is known as the free fermion limit of the XX Z chain, and it is the most studied 
one [25H28]. 

In the next section II we shall consider the Ising limit corresponding to the limit 
A — > — oo. The thermal correlation functions of the ferromagnetic string operator are 
obtained in Section III. Section IV presents several results for the four-vertex model which 
is the A — > — oo limit of the six-vertex model. In Section IV it is demonstrated that the 
effective Hamiltonian commutes with the transfer-matrix of the four-vertex model, the 
orthogonality of the eigen-functions is proved, and the relation between the four-vertex 
model and the boxed plane partitions [12J is discussed. The low temperature asymptotics 
of the correlation function is studied in the Section V. Discussion in Section VI concludes 
the paper. 
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2 The Ising limit of the XXZ chain 



Less studied limit of the XXZ model is the Strong Anisotropy (SA) limit A — > — oo 
jU|29H3T]. This limit is called the Ising limit since the simplest Hamiltonian in this limit 
turns out to be the Hamiltonian of the one-dimensional Ising model p[]: 

i _ _i M 

Jm^ - H xxz = H ls = — J>tn^ " 1) • (9) 

k=0 

However we shall study the strong anisotropy limit using the effective Hamiltonian Hsa, 
which is formally equivalent to the XX Hamiltonian flS}, provided the latter is supplied 
with the requirement forbidding two spin "down" states to occupy any pair of nearest- 
neighboring sites: 

^ M M 

Hsa = --J2 V & k+hk + ha z k )V, V = ]J(I - q k+1 q k ) . (10) 

fc=0 k=0 

The projectors V "cut out" the spin "down" states for any pair of nearest-neighboring 
sites. The local projectors onto the spin "up" and "down" states are of the form: 

q k = j ( E - a t) , Qk = \ (1+ Ofc) , & + ?k = I, k e M , (11) 

where I is the unit operator in (M + l)-fold tensor product of the linear spaces C. The 
point is that the Hamiltonians (Q and (fTPl) are commutative. Indeed, the Hamiltonian 
can be represented by means of (TTT|) as follows: 

M M 

H Is = N - Qk+iQk , N = Y,Qk, (12) 

k=0 k=0 

where N is the operator of number of particles commuting with the Hamiltonians Hxx 
and Hsa- On another hand, the definition of the projector V yields: 

M M 

V qk+iqk = ^ Qk+iqkV = 0. 

k=0 k=0 

Then, [Hsa, Hi s ] = and so the Hamiltonians ([9]) and ( JTUj) possess a common system of 
the eigen-functions. 

The wave function (j3J) takes the form in the limit A — > — oo: 

xf(u)=det(^- JV+fe) ) 1 <, fc < w J] M-O' 1 ' (13) 

l<n<l<N 

where the coordinates of the spin "down" states form a strict decreasing partition /x, as 
in ([3]). It follows from (TT3|) that the wave function is not equal to zero if and only if the 
elements /Xj, 1 < i < N, satisfy the exclusion condition: fa > fa + i + 1. In other words, in 
the considered limit occupation of nearest sites is forbidden, and the hard-core diameter 
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equal to duplicated inter-site separation arises. The exponential parametrization v3 = e ldj 
brings the Bethe equations to the form: 



N 



e i(M+l-N)8 k 



■if- 1 Y[e ld > l<k<N, (14) 



which is respected by the solutions 

W+l-N 



1 N 

where P = — 9j, and the Ij are integers or half- integers, depending on whether N is 

odd or even. It suffices to consider a set of N different numbers Ij satisfying the condition 
M — N > 1\ > J 2 > • ■ • > In > 0. The ground state of the model is defined by the 
solutions 

with the condition P = fulfilled. The notation 6 for iV-tuple (9i,02, ■ ■ ■ ,&n) of the 
solutions (1 15)) will be especially convenient further in order to stress that one is concerned 
with the solution of the Bethe equations. Otherwise, it is appropriate to use u as an 
indication that an arbitrary parametrization of the wave- function (TT3j) is meant. The 
eigen-energy of the ground state of the model calculated by means of (J2J) and ( 1T6|) is: 

7T IT N 

E N (0 Y ) = hN - cosec— — sin— — -. (17) 

v; M+l-N M + l- N K ' 

The eigen-energy of the Ising Hamiltonian ()9]) has the form: E\ s = N. 

The two limits, flS} and (fTUj) . are similar in the sense that their wave functions are 
expressible through the Schur functions [12] : 

q i \ q / x _ det(x Afc+jV ~ fc )i< J - fc < jV 

5 A (x) = Sx(x 1 ,x 2 ,...,x N ) = N _ k , (18) 

aet{Xj >i<j,k<N 

where A denotes the partition (Ai, A2, ■ ■ ■ , Ajy), which is an A^-tuple of nonincreasing 
non-negative integers: L > Xi > X 2 > ... > A^ > 0. Indeed, any strict partition 
M > fii > fi 2 > ••• >/ijv>0 and non-strict partition M + l — N > Xi > X 2 > 
■ ■ ■ > Xn > (denoted as \i and A, respectively) can be related by the formula Xj = 
fij — N + j, where 1 < j < N. In other terms, A = /x — 8, where 6 is the strict partition 
(N — 1,N — 2,...,1,0). Any strict partition \i with the elements respecting the exclusion 
condition /ij > + 1 is connected with the non-strict partition A by the relation 
A = n - 25, where M + 2(1 - N) > X 1 > X 2 > ■ ■ ■ > X N > 0. So, the wave function 
f|T3|) may be represented as X^, A ( U ) = ^(u 2 )? an d the corresponding state- vector ([3]) is 
specialized as follows: 

M 

IMu))= E 5 a(" 2 ) IlK") efc ft) ■ ( 19 ) 

AC{(Af-2(AT-l)) JV } fc =0 
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Summation in (TT9|) goes over all non-strict partitions A into at most N parts so that each 
is less than M - 2(N - 1), i.e., M - 2(N - 1) > Ai > A 2 > ■ ■ • > \ N > 0. Summation 
over the non-strict partitions A is equivalent to that over the strict partitions with the 
elements respecting the conditions /ij > + 1 and A = \i — 28. 



3 Survival probability of the ferromagnetic string 

We shall study the thermal correlation function of the states with no spins "down" on the 
last n + 1 sites of the lattice (would be called as survival probability of the ferromagnetic 
string) defined by the ratio: 

{*»(»•) I ^» j}}-n 

where Hq,a and Tl n are defined, respectively, by (TTU]) and (ITT]) , and (3 is inverse temperature. 
Besides, 6 V indicates that the eigen-states are calculated on the solution (TT6]) . Our calc- 
ulations will be based on the following relation for the Schur functions (TT8|) which is due 
to the Binet-Cauchy formula [32J: 

£ S x (xl,...,x 2 N )S x (yl,...,y' N ) = det ( i^g^ ) 

>< n fe 2 -^)" n w-a&r 1 - (21) 

l<k<j<N l<rn<l<N 

Summation in (1211) goes over all non-strict partitions A into at most N parts so that each 
is less than L: L > Xi > A 2 > ■ ■ • > > 0. The notation for the range of summation in 
( 12~T|) will be extensively used in the rest of the paper. 

Firstly, it is straightforward to calculate, using (12T|) . the following average of the 
projector 77 n (fTUj) over the state- vectors (|T9|) : 

(*^(v) | n n | ^( U )> = Yl %(v- 2 )%(u 2 ) 

\C{(M~2N-n+l) N } 

det ^4t^ (22) 



V(U 2 )V(V- 2 ) V h<j,k<N 



3 

where the notation for the Vandermonde determinant is used: 



V(u 2 ) ^ J] ^ - u "m) ■ (23) 



l<m<«<Af 



Right-hand side of f )22p . taken at n + 1 = 0, gives the answer for the scalar product of 
the state- vectors (I19p . i.e., for jv ( v ) | ^iv( u ))- Equation (12"2"|) . taken on the solutions 
fn6|) . allows one to obtain the so-called Emptiness Formation Probability E(n,N), which 
expresses the probability of formation of a string of the empty (i.e., spin "up") states on 
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last n + 1 sites of the lattice 
E{n,N) = 

+ 





tin 











\ _ e in{ki-kj) 



M + 



— det((l 



r?. 



A/ 



M-iV + 



7^ 



M+l-iV 



(M - JV+ 1)(1 - e*(*'-^)) v "V i<«,j<jv 

The solutions of the Bethe equations ( 1T3|) constitute a complete set of the eigen-states 
and therefore one can obtain the resolution of the identity operator: 

1 



I 



(M + 1)(M + 1 -N) 



N- 



rl>(< 



i0\\ 2 



(24) 



where the summation over {0} implies summation over all independent solutions (TT5l) . 
Here the exponential parametrization for the solutions of the Bethe equations is used in 
the compact form, e l9 = (e 101 , e 102 , . . . , e l0N ). We shall calculate the nominator of (120 j) 
using insertions of the resolution ( l24"j) . Taking into account that 



(*iv(v) | e"^ A |**(0)> = (\Mv) |^ JV (0))e-^W 5 

we obtain: 

(^(v) | n n e-^n n |*jv(u)) 



i 



(M+l)(M + l-iV) 7V - 



rE 



-PE N {0) 



{0} 



x |V(e ie )| 2 p(v- 2 ,e lW )P(e- iW ,u 



-2 



where 



(25) 



(26) 



The range of summation in ( 126]) is taken as follows: A L , A^ C {(M — 2N — n + 1)^}. 
Then, using the relation (j2ip we obtain from (I25p : 



x det 



(*jv(v) \n n e-e H ^n n \* N (u)) = 
1 

(M + 1)(M + 1 - iV)^ 1 V(u 2 )V(v- 2 ) 
1 



E 



-PE N {6) 



f^ e ie ivT 2^M-N-n+l 

1 - e**uj" a 



det 



M-AT>/i>/ 2 ->i'jv>0 
1 _ [ u 2 e -i9 P \M-N-n+l 



(27) 



±<ij<N 



1 — ufe~ i0 p 



l<l,p<N 



where summation goes over the ordered sets {Ik}i<k<N that parametrize (TT5"j) . and 
Ejsf(0) is given by (|17|) . Expression for T(6 v ,n, (5) f l2U|) appears from f[2"T|) as follows: 



T(0 v ,n,/3) 



1 



(M+l) 2 (M + l-iV)^- 2 
x det 



E « 

M-Ar>/i>/ 2 ->/]v>0 
1 _ ^(M-TV-n+l)^-^) 

1 _ e i(fi-^) 



P(E N {0)-E N {O v )) 



l<p,l<N 



(28) 



where # v implies the ground-state solution (TTB"]) . and En{O v ) is the ground-state energy 
in the limit of strong anisotropy. 
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4 Four-vertex model and boxed plane partitions 



4.1 The XXZ Hamiltonian at A 



-oo and the four- vertex model 



The six-vertex model on a square lattice is defined by six different configurations of arrows 
pointed both in and out of each lattice site. A statistical weight Wk (k = 1,2, ... ,6) is 
ascribed to each admissible type of the vertices (Fig. 1). Representing the arrows pointing 
up or to the right by the solid lines one can get the alternative description of the vertices 
in terms of lines floating through the lattice sites. So far the bonds of a lattice may be 
only in two states - with a line either without it, there is a one to one correspondence 
between the admissible configuration of arrows on a lattice and the network of lines - the 
nest of lattice paths. The four-vertex model is obtained in the limit when the weights 
W\ = W3 = 0. 






Wi 



Vt'2 



M'3 



VV'4 



Wi 



MY, 



++++++ 

Figure 1: The vertex configurations of the 6- vertex model. 



The L-operator of the six- vertex model is equal to |36| : 

L 6v (n\u) = 



we 7 °™ — u 1 e 7 °™ 2 sinh(27) <y n 
2sinh(27)a+ ue^< + u' 1 e^< J' 



where h6C and 7 = |ArchA, and it satisfies the intertwining relation: 

R(u,v)(L &w (n\u) (g> L ev (n\v)) = (L 6v (n\v) (g> L 6v (n\u))R(u,v) 
Here, R(u, v) is the (4 x 4)-matrix: 



R(u, v) 



( f(v,u) \ 

g{v,u) 1 

1 g{v,u) 

V f(v,u) / 



where 



f(v,u) 



V? — V 2 



g[v,u) 



uv 



V? — V 2 



( e 27 _ e -27) ; 



and u, v e C. 



Let us consider the following transformation of the L-operator 
L(n\u) = e h <e {uj/2)aZ L 6v (n\u)e- {uj/2)aZ 



u e v 



2 sinh(27) e u + h < o7 



2 sinh(27) e ~ w+h < cr+ ue ( h -~f)< + u -i e (MiX 
8 



(29) 



(30) 



(31) 



(32) 



(33) 



The operator (133|) is intertwined by the transformed .R-matrix of the form: 



R(u,v) = (l®e- haZ )R(UiV)(l 

( f(v,u) 

g(v,u) 



e h °°) 



(34) 






V o 





2h 



e 2h g(v,u) 







\ 



f(v,u) ) 



If one puts h = u = 7, then the limit 

lim e~ 21 L(n\u) = L 4y (n\u) 



7—^00 



-uq n a' 



U 1 q n 



(35) 



where q n = <y\o~ (see fTTT]) ). gives us the L-operator of, so-called, jour-vertex model. 
Analogously, the limit 7—7-00 transforms the matrix e~ 2j R(u, v) into the .R-matrix of the 
four- vertex model |22|: 



R(u, v) 



( f(v,u) \ 

g(v,u) 1 

g{v,u) 

V f{v,u) J 



(36) 



where 



f(v,u) 



u 



u 2 — V 2 



g(v,u) 



uv 



u 2 — V 2 



(37) 



The monodromy matrix of the models is defined as the matrix product of L-operators: 
T{u) = L(M\u)L(M - l\u) ■ ■ ■ L(0\u) = I ^ ^\ ) . (38) 



The transfer matrix is the matrix trace of the monodromy matrix: t(u) = TrT(w). 
In |33| it was proved that the transfer-matrix of the six-vertex model commutes with the 
Hamiltonian of the XX Z model ([T]): [H 6v , t 6v (u)} = 0. 

The transfer-matrix of the four- vertex model r 4v (-u) satisfies the property: 



Vt 4v (u) = r iv (u)V = t 4v (m) 



(39) 



where V is the projector defined in (fTOi) . Indeed, using the explicit expression for the 
entries of the product of two //-operators, 



L(n + l\u)L{n\u) = ( ^ n+1 ' n ^ n+1 > n ) 



(40) 



one can show that 

C n+ i tn L{n + l\u)L{n\u) = L{n + l\u)L{n\u) C n+ i^ n = L{n + l\u)L{n\u) , (41) 

where C n+1 , n = 1 - q n+ iq n . It follows from flU]) that C n+ i, n T 4v (u) = t 4v (u) C n+ i, n = t 4v (u) 
at n = 0, . . . , M — 1. Invar iance of the matrix trace with respect of cyclic permutation 
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implies that Co.a/T4v(w) = T4 V (k) Co,m = r 4v(^), and therefore the relation fl39|) is also 
fulfilled. 

Further, one can advance the following 

Proposition: The transfer matrix of the four-vertex model commutes with the XX Z 
Hamiltonian taken in the Ising limit: 

[#SA,T 4 »] =0. (42) 



Proof: In order to prove (l4z|) . it is suffices to consider Hsa (|T0|) at h — 0. We shall 
use the method proposed in [33] and then developed in [3l]. Generalizing this approach 
with respect of the problem under consideration, we shall take into account the relation 
(subscripts are omitted): 

C n +i, n \h n+1>n , L(n + l\u)L(n\u)] C n+ i <n 

(43) 

= C n+ltn (Q(n + l\u)L(n\u) - L(n + l\u)Q(n\u))C n+ i, n , 

where the operators h n+ i, n and C n+ i. n are introduced in ([8]) and (HIT) , respectively. If so, 
we obtain the vanishing of the commutator in question, 

- AI . A/ 

n=0 n=0 

provided that 

p n+1 , n ,T(u)] = L{M\u) ■■■L(n + 2\u) {h n+ltU ,L(n + l\u)L{n\u)] L(n - l\u) ■ ■ ■ L(0\u) . 

(45) 

Furthermore, the commutators of the entries from the right-hand side of ( 14U|) look as 
follows: 

[^n+l,n) ^4n+l,n] = ~ [^-n+l,n) ^n+l,n] = Qn+1 ~ Qn 

\h n+ltTl , B n+lj7l ] = -ua~ +1 + u~ x a~ + ua~ +1 q n - u~ l q n+1 a~, (46) 
\h n+ x jn , C n+1)n ] = -tt _1 cr+ +1 + Mff+ + u -1 cr+ +1 g„ - uq n+1 a+ . 



In turn, one obtains from (146]) : 

C"n+l,n [^-n+l,n) ^n+^n] C n j r \^ l C n +i j?1 (<Jn+l ^nj^n+l,?!; 

C-n+l,n \hn+l,m -Sn+l,n] C n +i jn Cn+l,n 

C n +l,n [^-n+l.n; Cn+l,n] C n+ i jn = U 1 0" r [ +1 + ' ucr ^) Cn+l,m 

C"n+l,n [^-n+1,71) ^n+l,n] ^n+l,n ^-'n+l,n( *?n+l "I - '?7i)^'Ti+l,n- 

On the other hand, using the diagonal matrix 

Q(n \u) = 
we obtain the following relation for the L-operators fl35|) : 



(47) 



()(/, + i|„ )!(/,(, )-i(„ + i|»K?(„|») = ( J?+l l n Ua 7 1+ ^ )-( J8 ' 

~~ " °n+l + ucr n — Hn+1 i (/n 
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Comparison of ()47|) and (l4~8j) demonstrates that the relation ( l4~3j) . as well as the statement 
of the relation ( 14"2|) are valid indeed. ■ 

The state vector of the four-vertex model is constructed in the framework of the 
algebraic Bethe ansatz as follows: 



N 



\y N (u)) = l[B Av (u i )\t}. (49) 



i=l 



This vector is the eigen-vector both of the transfer-matrix of the four-vertex model and 
of the Hamiltonian ( II Up . provided that the parameters Ui fulfill the Bethe equations 0141) . 
It was shown in [22J that this vector can be represented in the form ( fl~9]) . 

According to the representation (122]) taken at n + 1 = 0, the scalar product 
(^^(v) | \1/Ar(u)) is proportional to the determinant of the matrix, say, {Tkj)i<k,j<N- 
The Bethe equations (lT4"j) enable one to represent this determinant as follows: 



det(T fei ) i< fcJ <jv = det ( ^ Z^TTp. ) = ( II ) det (^fcj)i<fc,j<iv , (50) 



2\A4-2V+2\ / N x N 



J / \fc=l 

where the entries Akj take the form: 



■ M+3-V 



N N 



- 7 ? n «?- n »n- («) 

"fc U J V=l,i?4fc 1=1, lytj / 

It is crucial that the determinant of the matrix {Akj)i<kj<N vanishes provided the par- 
ameters u and v are independent solutions of the Bethe equations 0T4j) . Really, this 
matrix has a non-trivial eigen-vector \j with zero eigen-value: 

N N 2 2 

E^ = °' ^ = K 2 -^)n^— i- ( 52 ) 

3=1 i=i,¥i U i Ul 

Validity of (152|) is justified by two identities: 

N / V 

\ ^ / yfc - a^fc 1 r Uk - xi 

N / AT \ N 

Therefore, the scalar product (^^(v) ^^(u)) vanishes (i.e., the state-vectors are ortho- 
gonal) provided the sets of the parameters u and v represent independent Bethe solutions. 
Completeness of the system of the state- vectors was proved in [HI5U]. 



4.2 Boxed plane partitions 

Consider the four- vertex model [22]. There is one to one correspondence between the 
strict plane partitions (i.e., the plane partitions that decay along each column and each 
raw [12J) and the nests of the admissible paths on a square lattice of the size 2N x (M + 1) 



11 




Figure 2: The admissible nest of lattice paths of the four-vertex model and the corre- 
sponding strict plane partition. 

with the following boundary conditions: all arrows on the left and right boundaries are 
pointing to the left, while the arrows on the top and bottom of the first iV columns 
(counting from the left) are pointing inwards and the arrows on the top and bottom of 
the last iV ones are pointing outwards (see Fig. 2). 

It was proved in [22J that if one puts Vj = q~i and Uj = q~^~ in the scalar product of 
the state- vectors (H9l). then 



,g~f) | ¥*(!,..., g^)) = q~^^ Z? p (iV, N, M) 



(53) 



where Z^ PP (N,N,M) is the generating function of strict plane partitions placed into a 
box of size N x N x M: 



Z° PP (N, N, M) = q 



_ „iV 2 (iV-l) 



l<j,k<N 



I _ gM+3-j-k 



- qj+k-i 



(54) 



Being taken at q = 1, this formula gives the number of strict plane partitions in a box of 
size N x N x M: 



Z^N, N, M) 



It is straightforward to obtain 



n 



M + 3 - j - k 
3 + k-l 



l<j,k<N 

that the entry 



lim (^jvfg 



n n |^(l,...,g— )) = ffi(JV,JV,M-n-l) 



(55) 



(56) 



is equal to the number strict partitions in a N x N x (M — n — 1) box. 



5 Low temperature limit 



Let us obtain the low temperature estimate of the correlation function (|2T)|) . The survival 
probability of the ferromagnetic string is written as follows: 

|2 



T{0\n,P) 



\V(e l 



V 2 



J2 e -P(E N (e)-E N (ev)) | v(e<9) v{e -i9^ e 

{0} 



J0 V \\ 2 



(57) 
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where V = (M+ 1)(M + 1 — N) N 1 . Here V(e %e . e l ° v ) is given according to fl26|) . provided 
the summation domain is given by A C {(M — 2N — n + 1)^}. 

If the chain is long enough while the number of quasi-particles is moderate, i.e., iV <C 
M, we replace the sums in (1571) by the integrals as follows: 



r ^ V ^^^(^Ww II 2(l-cos(^-^)) 

^ ' l<r<s<N 

x fi(r d A)/E^-™°» lne -» >e ^ n | e *,_ e «f. (58) 

i=l ,/u l<k<l<N 

In the considered limit, we may asymptotically put cos 6^ ~ 1 and cos6j w 1, V7. At the 
large enough /3, we approximate (158]) in the following way: 

TV 2 ('\ 1") /?\N 2 /2 _ 

^"^^ wQ n < r - s ) 2 . ( 59 > 

V ; M l<r<s<7V 

oo oo oo jV 

I = / ■ ■ ■ e 1=1 ' J J Jxfc — xj| dxidx 2 ■ ■ ■ dxjsr . (60) 

~ ~ "L Kk<l<N 
— oo— oo — oo — — 

The integral X (160]) is the Mehta integral, and its value is known [35J: 

7r Af / 2 r(i|iv + 1) 

~~ 2N(N-l)/2 ' 

where 

N 

II |r- s | = r(i|jv), r(i|N) = JJr(n). 



l<r<s<iV n=l 



Eventually, taking into account (156]) and (|59[) . we express the answer in the following 
form: 

r SA (0,n^) * (Z^(iV,iV,M-n-l)) 2 ^, Cjv _ 

The low temperature decay of the survival probability of the ferromagnetic string is gov- 
erned by the critical exponent N 2 /2, while its amplitude is proportional to the squared 
number of strict plane partitions in a box of size N x N x (M — n — 1). 



6 Conclusion 

In this paper we studied the XXZ Heisenberg magnetic chain in a specific limit of the 
anisotropy parameter: A — > — oo. The, so-called, survival probability of the ferromagnetic 
string has been calculated over the iV-particle ground state of the model. It was proved 
that the Hamiltonian of the Ising limit, previously introduced in the papers [21 ] 129 ] 130]. 
is correct as the limit A — > — oo of the XXZ Hamiltonian. Finally, we have found 
the leading term of the low temperature asymptotics of the survival probability of the 
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ferromagnetic string for the model in the finite volume and at fixed number of quasi- 
particles. We have demonstrated that the amplitude of the asymptotics is proportional 
to the squared number of strict plane partitions in a box, while its critical exponent is 
proportional to the squared number of particles. Comparing this result with that obtained 
in [SUET], one can conclude that this type of the behavior is universal for a special class 
of the integrable models. 
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